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We study the behavior of magnetization curve as a function of magnetic field in the immediate 
vicinity of the magnetization plateaus of ID electron systems within the bosonization formalism. 
First we discuss the plateau that is formed at the saturation magnetization of ID electron system. 
Interactions between electrons we treat in the lowest order of perturbation. We show that for isolated 
systems, where total number of electrons is not allowed to vary, magnetic susceptibility stays always 
finite away of half filling. Similar statement holds for many other magnetization plateaus supporting 
nonmagnetic gapless excitations encountered in ID electron/spin systems in the absence of special 
symmetries or features responsible for the mode decoupling. We demonstrate it on example of the 
plateaus at irrational values of magnetization in doped modulated Hubbard chains. Finally we 
discuss the connection between the weak coupling description of saturation magnetization plateau 
and strong coupling description of zero magnetization plateau of attractively interacting electrons/ 
antiferromagnetically interacting spin 1 Bosons. 



Magnetization process of ID electron systems is fairly well understood theoretically. For purely spin systems (where 
charge fluctuations are frozen) the magnetization curve (magnetization plotted versus magnetic field) can show plateau 
like behavior (at zero temperature) for certain rational values of magnetization^ where magnetic excitations develop 
gap^. 

Experimentally different materials have been sinthesized exhibiting plateaus in their magnetization curve, and 
having magnetic structure believed to be modeled by quasi- one dimensional geometries. For example dimcrizcd spin 
5 = 1 chain^ exhibits magnetization plateau at 1/2 of saturation magnetization as predicted by theory^. Other 
quasi-lD compounds showing plateau at half of the saturation magnetization are spin S ~ 1/2 compounds^i^. 

Phase transitions at the edges of the magnetization plateau is usually described within the commensurate- in- 
commensurate universality clasa^ii. As a hallmark feature of commensurate - incommensurate phase transition the 
magnetization shows square root dependence on magnetic field in the vicinity to the edges of the plateau, giving rise 
to infinite magnetic susceptibility at the edge points. The simplest examples include: magnetization plateau of spin 
S chain at the saturation magnetization and a plateau of integer spin 5* antiferromagnetic Heisenberg chain (or even 
leg spin S ladder) at the zero magnetization and we depict them on Fig.(l). 

When doped with charge carriers the plateaus can appear also at irrational (doping controlled) values of 
magnctizatior*^. Interesting property of such plateaus is that they can support nonmagnetic gapless excitations. Other 
examples of plateaus with similar feature may be found in frustrated spin chains, like trimerized S = 1/2 zig-zag chain 
at 1 /3 of saturation magnetization where gapless chiral excitations coexist with gapped magnetic excitations^, or spin 
S generalization of the integrable t — J chain doped with S — 1/2 carriersi^. 

In this work we investigate in detail the edge points of plateaus with gapless nonmagnetic excitation where modes 
describing magnetic (gapped at the plateau) and non-magnetic (gapless) excitations do not decouple. We will argue 
that in such situations, in drastic contrast to the square root behavior characterizing the commensurate- incommcn- 
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FIG. 1: magnetization curve of integer spin S Heisenberg chain. 
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FIG. 2: Portion of the magnetization curve of a doped dimerized Hubbard chain. 



surate phase transitions, magnetization depends linearly on magnetic field at the edges of plateau, thus rendering 
magnetic susceptibility finite at the edge points. 

The paper is organized in the following way: First we describe within the effective theory bosonization approac h^ 
the simplest plateau of saturation magnetization encountered in free electron system. Next we add interactions 
between electrons and treat them in the first order of coupling constants. We compare our results to exact results 
available for the Hubbard model for both positive and negative Hubbard couplings. In the case of attractive Hubbard 
interaction we corrected the missprint in the exact expression of the saturation susceptibilit y^'^i and convinced that 
bosonization recovers exact expression of saturation magnetization in the lowest order of onsite Hubbard interaction 
for both signs. We apply the similar theory to other magnetization plateaus that are encountered in more complicated 
problems, supporting gapless excitations, and in particular to the plateau in dimerized doped Hubbard chain at doping 
S = 1 — n, where n is a filling. As was shown in Refi^ at magnetization m = S/2 plateau is formed and we show 
that magnetization changes linearly with the field once plateau is closed by magnetic field on both sides as depicted 
on Fig. (2). Finally we discuss the connection between the weak coupling description of saturation magnetization 
plateau and strong coupling description of zero magnetization plateau of attractively interacting electrons. 

In appendix we give our bosonization rules and for completeness give exact expressions for saturation susceptibility 
for Hubbard model for both values of Hubbard coupling constant. 



II. MODEL AND EFFECTIVE FIELD THEORY 



The microscopic model, general enough for our purposes is described by the following lattice Hamiltonian of the 
one dimensional electrons: 

N N 
N ^ N 

+V^^ni^ni+i„ - -^(?i^,T -rii.i) (1) 

i — l (7 2—1 

We will always work at the fixed filling. Above 7ii_o- = c|^Ci^cr is the number operator of spin a ~ {t, i} electrons at 
the i-th site, U is the onsite coupling constant, V describes interaction between the same spin electrons on neighboring 
sites, h is an external magnetic field, ti ~t + 1)' is modulated hopping amplitude with period 2 (dimerized). 

We will study the continuum limit of the above model using bosonization approacb^^*^. Bosonization rules are 
summarized in Appendix A. 

In bosonic formulation effective field theory of free fermions (putting t' = U = V = Q in Eq. ([T])) is just a direct 
sum of two decoupled Gaussian models, one for each spin component: 

^=-Ey/ dx[{d,<^,f + {dj,f]. (2) 

where velocities are determined by the fillings: Va = 2tasanma, with a standing for lattice constant that wc will put 
to 1. 
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III. SPIN SUSCEPTIBILITY AT THE SATURATION MAGNETIZATION 



Magnetic susceptibility of the Hubbard model with fixed filling is always finite at the saturation magnetization for 
any interaction strengtbi^iiiii^ (including U > 0, U < 0, and U = 0) except for half fillingi^ and empty or completely 
filled bands. It is instructive first to develop bosonization analysis and reproduce the above result within the effective 
field theory language. This will serve to us as a warm up procedure for more complicated cases considered later. In 
this section we put <' = and t = 1 in Eq. 



A. Bosonization for free electrons at saturation magnetization 

First we put U = = in Eq. ([T]) and apply bosonization for free electrons at saturation magnetization. At 
saturation magnetization band of down- spin electrons gets empty, while the band of up- spin electrons gets either 
partially filled, if filling is less than half, or completely filled exactly at half filling. To mimic the bottom of the down 
spin electron band we add an appropriate relevant sine- Gordon potential to free bosonic theory describing the down- 
spin electrons (see Appendix B): 

^ = Y [idtd^ir + id.^lf] + y [(5*0t)' + (d.H)'] 

1, . ,dr(t>\ r- Ah dT(t>r 

+ - {hw + Ah) ^^ + W cos VIttcPi-— — ^ 

2 y/TT 2 y/7T 

- A,(^ + ^) (3) 

Above = 2sin(7rn), is the velocity of up spin band at saturation, hw closes the gap opened by the cosine term 
and Ah is a deviation of magnetic field from the saturation value. Down- spin electron parameters and W are 
selected in such a way to mimic optimally the dispersion at the bottom of the down spin electron band. As it will 
become obvious we will not need their explicit values for calculating the magnetic susceptibility at the saturation 
magnetization. The chemical potential is a function of magnetic field (Lagrange multiplier enforcing constraint) and 
keeps the number of electrons fixed. 

In bosonization the constraint of fixed particle number looks as: 

(a,^4+9,</)T) =0 (4) 

On the other hand uniform magnetization density is given by the following expression: 

Am = 'f-'^^^ (5) 

2V7r 



At the saturation magnetization wc have from Eq. ([5]): 

(d^cj)^ ) _ Ah/2 + An{h) _ A/i/2 + A/z(/i) 



(6) 



y/n nvi 2tt sin im 

Using Eqs. ^ and ([5]) from Eq. ^ we can express the Lagrange multiplier A^ as follows: 

A^i^TTv^Am- Ah/2 (7) 

In addition using the sine- Gordon Hamiltonian for down- spin electrons we have second equation^ii (for details see 
Appendix B): 



An, = Am^ -^/\Ah/2^^A^J\ 



f 



TT 



-J\Ah/2- {T:v^Am- Ah/2)\ (8) 



where for the last equality we have used Eq. ([7]). The correct root of the Eq. ([5]), one vanishing with A/i, is: 

Am=^-(^-fO(A/.3) (9) 

Yielding for the magnetic susceptibility at the saturation magnetization: 

X=l^\f^sat= = 7, • (10) 

ah 7rw| 27r sm Trn 
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B. Including interactions perturbatively 

From the Eq. (fTO|) we see that the susceptibihty is finite except at half filhng. Important observation is that 
this result is robust to interactions, one can treat interactions between electrons perturbatively, because the saturation 
magnetization plateau is caused by the band effect and exists already at the free electron level. To demonstrate this we 
will include perturbatively Hubbard coupling U and next nearest neighbor interaction V and obtain exact expression 
of susceptibility in the first order of these coupling constants. Interactions U and V in bosonization formulation close 
to saturation magnetization are translated as: 

U V V 

—d.^(l)ld^(j)^ + —d^(j)id^(f)i + —d^4>^d.^4>^ (11) 

TT TT TT 

Other terms are either incommensurate- meaning that they rapidly oscillate and cancel out in the continuum limit or 
irrelevant (like Umklapp process at half filling for the attractive case) and thus are neglected. 

Using constraint ^ from ([5]) one can express magnetization by d^ij)^ simply by minimizing the energy: 

Where H is given by the sum of ([3]) and (fTTj) . This fixes magnetization in the ground state configuration as follows: 

Am = ■^=— = (13) 

Ato(1 ) = (14) 

= 7rm Ato(1 - ^ ~ ^^ ) - Ah/2 (15) 

Now we need another equation to exclude Afi. We adopt a mean field decomposition of interaction terms and for the 
second equation we obtain: 



Am -^/\Ah/2 -An+{U- 2V)Am\ 

TT 

= -J\Ah-Am{nv^ - 2{U - 2V))\ (16) 

TT » 



Am = ^ + 0{Ah'^) (17) 



where 



For the susceptibility at the saturation magnetization in the first order of coupling constants we obtain: 



Z7rsm(7rri,) \ TrsmTm/ 

The plus and the minus signs in front of U and V respectively are intuitively clear: repulsive Hubbard interaction 
facilitates to fully polarize the system while repulsion between the similar spin components hinders it. Moreover 
number of nearest neighbors along the chain is twice the number of sites, hence contribution from V is multiplied by 
factor of 2 in comparison with U . 

In the case of the pure Hubbard model {V ~ 0) we corrected a misprint of factor of 2 ii^i^ (one can repeat Bethe 
Ansatz calculation close to saturation magnetization and get convinced that bosonization indeed gives correct result 
in the first order of U for both signs, see Appendix C). 

The above effective theory developed at the edge of the saturation plateau is not only valid in this particular case, 
but also applies to other plateaus, e.g. those that arc characterized by gapless non-magnetic excitations which we will 
discuss in the following section. 
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IV. PLATEAU OF DOPED DIMERIZED HUBBARD CHAIN AT IRRATIONAL VALUE OF 

MAGNETIZATION 

In this section we will apply our bosonization description that we developed in previous sections to calculate magnetic 
susceptibility at the edges of plateau of doped dimerized Hubbard chain {V ~ 0,t' ^ 0) at doping dependent value 
of magnetizuation 2m = S^. This plateau shares the similar physics of the saturation magnetization one. The exact 
equivalence with the saturation plateau is due to the fact that description is perturbative in interaction parameter, 
plateau exists already at the level of quadratic in Fermi fields Hamiltonian (due to the modulation of exchanges) and 
the limit U ^ is non- singular. 

The bosonized Hamiltonian at the plateau of dimerized doped Hubbard chain, (first setting Hubbard interaction 
to zero) is given by Eq. ([3]), since at this plateau one degree of freedom (up- spin) remains gapless, whereas down- 
spin electrons meet the commensurability condition and thus are gapped. Now we can add perturbatively Hubbard 
interaction, or more general interactions, and follow the steps that we developed in previous section. In this way we 
will obtain susceptibilities in the lowest order in interaction constants. Situation is manifestly similar to the saturation 
magnetization case. From this equivalence we can answer a question^ on the critical properties at the edges of these 
plateaus. Namely, magnetization will increase always linearly with the magnetic field, unless the total number of 
electrons is allowed to vary. If no constraint on the total number of particle is imposed, plateau itself will disappear—. 
Once again, due to the linearly dispersive degree of freedom, potential Van Hove singularity from the bottom of the 
gapped mode is suppressed by the constraint of keeping total number of electrons fixed. 

We think that similar effective theory can apply to the doping dependent magnetization plateau in an integrable 
spin -S generalization of the t — J chain doped with spin 5—1/2 carriersi^. As explained within the Bcthe Ansatz 
methocU^ in this model quantum numbers of excitations carrying spin and charge get coupled at the right edge of the 
doping dependent plateau and this can cause the linear increase of magnetization as opposed to the naively anticipated 
square root behavior. 

V. SPIN SUSCEPTIBILITY AT THE ONSET OF MAGNETIZATION 

More subtle situation arises when magnetization plateau is formed due to the interactions (four fermion involving 
terms), and to describe correctly the edge of the plateau one needs to use non- perturbative approach. One example 
of such situation is the weak coupling limit of Hubbard model of attractively interacting electrons where the effects 
of curvature of the Fermi surface on commensurate- incommensurate transition in spin sector away of half filling were 
emphasizedii. However in the dilute limit of the same model a simple approach was developedii using the reasoning 
along the lines suggested in^ which allows to obtain exact expression for susceptibility at the onset of magnetization. 

A. Attractively interacting electrons 

In this section we remind the bosonization description of zero magnetization plateau of attractively interacting two 
component fermionsii. In the dilute limit the similar approach as that at the saturation magnetization was developed, 
treating bound pairs and up spin electrons as independent particles, but constrained by the total number of particles 
being fixed. The independence could be justified by a dilute limit, and treating bound pairs as hard core bosons can 
be justified by the inherent strong coupling nature of the dilute limit. Thus in the dilute limit the heuristic model 
developed in strong couplingi^ is expected to yield exact expression for susceptibility. For completeness wc remind 
the calculations developed in Ref.[13l in dilute limit. The Hamiltonian is given by the following expression: 

{hw + Ah) 9^01 r- 

- ^ + W cosV 47r0| 

2 y/n 

- f^{—^ + —^+n) (20) 

Where index p stands for the bound pairs. The above Hamiltonian is manifestly similar to the one describing the 
effective theory at the saturation plateau except that magnetic field does not couple to the pairs. In the same way 
as outlined before one can easily work out the magnetic susceptibility. The dispersion of uncompensated up spin 
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FIG. 3: Variation of up- spin electron number in two interpretations: a) up and down spin particles, and b) up spin particles 
and pairs. Both parts correspond to the same configuration of spins, but due to different interpretations up spin electron 
variations differ by factor of 2. 



electrons at the critical field looks: 

i?T(fc) = y^^^fc^TA^-Ac.|^ (21) 

where A is a spin gap in the dilute limit, and U| is the velocity of uncompensated spins in the strong coupling (for 
Hubbard model, however their values do not appear in the expression of susceptibility. On the other hand bound 
pairs disperse linearly with velocity Vp ^ Q . Repeating analogous to the saturation case calculations we obtain for 
magnetization: 



Am = — \ —{h — her — VpirAm) (22) 

and for susceptibility: 

/C\n— const — ■ (^*^) 

TTVp 

For Hubbard model Vp ~ is the charge (soft) mode velocity, so we get: x\n=const = ^^r^- Thi is the low density 
limit of the exact strong coupling result: x\n=const = 27r^J(iL„)2 — ■ 

Note in the picture of up- spin and boimd pair particles Am = ^Arij = — A?ip here, different from the case where 
we considered the up and down -spin electrons where we used Am = A?i|. The differing 1/2 factor comes from a 
simple consideration sketched on Fig. (1) where one can see, that AA^|.b = 2AAf|.a- 

This approach captures the intuitive physics at the commensurate- incommensurate phase transition: the spin 
susceptibility is diverging because of the Van Hove singularity in the density of states of the uncompensated up- spin 
electrons when total number of electrons is not fixed. However, the linear dispersion of bound pairs, which also holds 
at the critical field, quenches the above-mentioned divergence when the total number of electrons is kept fixedi^. 
Certainly the above simple model provides with the correct physical picture, but it is valid in dilute limit, so it does 
not capture the special role of the half filled band. For the weak coupling limit close to half filling we refer to Ref . [13] ■ 



VI. SUMMARY 



We worked out effective field theory for the one dimensional electron systems describing the vicinity of the mag- 
netization plateau. First we treated free electrons at saturation magnetization and then we added interactions in 
the lowest order of coupling constant. We showed that away of half filling (for bipartite models that we considered) 
magnetization is a linear function of magnetic field at the edge of the saturation magnetization plateau if total number 
of electrons is not allowed to vary. The same theory describes also other than saturation plateaus, e.g. we described 
one encountered in dimerized doped Hubbard chain. Similar reasoning is easily generalizable to other plateaus sup- 
porting gapless nonmagnetic excitations in polimcrizcd doped Hubbard chains^. All these plateaus share the similar 
feature: they are adiabatically connected with non- interacting cases, that's why they are relatively easy to treat. 
More subtle situation arises when plateaus are opened due to the interaction, thus one can not describe them per- 
turbatively. For such situations we explained in details the connection between the zero magnetization plateau of 
attractively interacting two component Fermi system in dilute limifl^ and saturation magnetization plateau of free 
electron system. 
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The similar behavior of magnetization at the onset holds for spin 1 bosonsi^, where the spin degrees of freedom for 
zero magnetization are described by the 0(3) non linear Sigma model with gapped spectrmii^S, while charge degrees 
of freedom are gapless. Calculation simplifies in the dilute limit, which inherently falls under the strong coupling, 
where the system is made of strongly bound singlet pairs2ii2^ and the similar reasoning to that developed in Section V 
shows that magnetization increases linearly as soon as the magnetic field suppresses the spin gap if the total number 
of bosons is not allowed to vary. However to calculate for general situation magnetic susceptibility at the edge of zero 
magnetization plateau one should go beyond the spin- charge decoupling approximation (which follows from linearized 
hydrodynamics) in the effective formulation. 

Very generally when gapless excitations are present at the plateau and they are not decoupled from the gapped 
magnetic excitations due to special microscopic symmetrical^, magnetization of isolated system shows linear growth 
with magnetic field instead of the square root behavior. 

Our findings may be relevant for real quasi one-dimcnsional spin gapped materials, where one can observe experi- 
mentally crossover from square root dependence, at half filling, to linear growth of magnetization with magnetic field, 
away of half filling, in the vicinity of zero magnetization plateau. 
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APPENDIX A: BOSONIZATION CONVENTIONS 

Our bosonization rules are based on the following correspondence between the fermionic and bosonic operators^^: 



Where kp^„ — tm^ arc Fermi momenta expressed in terms of fillings: — < X^iLi "-i,'^ "Ha are Klein factors, 
ensuring fermionic commutation rules between electrons of different spin orientation: 

{Tla.lJa'} ^"iSaM' (A2) 

and and '^L,a are right and left chiral bosonic fields of spin a electrons with the following equal time commu- 
tation relations: 

1 

[*i?,a(a;), <^R,a'{y)] = 6a,cr'^sign{x - y) 

[^L,a{x),^L,a'{y)] = -S^,„, ^sign{x - y) (A3) 
From these chiral fields we construct the usual bosonic fields, and their dual counterparts: 

0aix) = $L,.(x) - $ii,<,(a;) (A4) 
with non- local duality relation between them: 

[Mx),Oaiy)] = te{y-x) (A5) 
where Q{x) is the Heaviside step function. Electron number operators are expressed in bosonization as follows: 

rii^a -* Ua + -^dx4>a{x) + — sm{2kpx + V^4>a{x)) (A6) 

v/tT TT 
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APPENDIX B: FREE ELECTRONS AT SATURATION MAGNETIZATION 

Here we give a brief review of the magnetization process of free electrons in the vicinity of saturation magnetization. 
For simphcity first we consider the grand canonical system, where the total number of electrons is not fixed. For 
external magnetic fields exceeding the saturation value, down- spin electrons band gets completely empty in the 
ground state. With decreasing the magnetic field from its saturation value the magnetization will also decrease. 
Change of magnetization will have two contributions: 1) from the partially filled up- spin electrons and 2) from the 
dilute down- spin electrons. The contribution from partially field up- spin electrons is linear with the field and it is 
given by Eq. ^ (where for grand canonical systems one should omit the Lagrange multiplier). Let us discuss the 
contribution from the down- spin electrons. Close to the saturation the density of down spin electrons is vanishingly 
small. In the absence of interactions the ground state energy density, given by its kinetic part is; 

eoi--^-^, (Bl) 

where ti^ = ^l/ L = An^ and mass was normalized to 1/2 (we remind that h ^ 1 throughout the paper). By definition 
chemical potential is given by the following expression: /i^ = 9eo|/9?i| which using Eq. (jBl[) gives: 



(B2) 



and accordingly: 



-1 



= ( ) = (B3) 



in particular we get that i^-l ^ ^■ 



For free fermions fi{T = 0) = Ep, N = L J^'' xd^, thus x is nothing else than density of states and its divergence at 
the bottom of the band is termed the Van- Hove singularity in the density of states. Most importantly this divergence 
is a universal feature - this result is not restricted only to free fermions: interactions (short range, falling quicker than 
l/x^) do not modify the above behavior. This one can check e.g. from exact solution of XX Z spin chain. Indeed 
in XXZ model, close to saturation, magnetization depends as a square root on deviation of magnetic field from the 
saturation value and numerical prefactor in front of the square root is independent of the value of magnetic anisotropy 
(Z part of the exchange). 

The relation between the particle density and chemical potential, given by Eq. (jB2p . holds also for massive 
relativistic fermions (Luther- Emery model), when chemical potential slightly exceeds the spectral gap, since the only 
ingredient needed for the above relation is the low energy dispersion, which must be quadratic in momenta - Eq. (|2ip . 
The Luther- Emery model on the other hand is equivalent to the sine-Gordon model, that is a model of real scalar 
field selfinteracting with cos/3</) term for particular value of /3 = ^/4^2L1^. This explains origin of Eq. ([8]) which 
follows from Eq. (|B2p . where for chemical potential one has to put: Ah/2 — Ajj., which couples to the density of 
down spin electrons given in bosonization by: dx4>i/\/TT. 

Now we discuss how the situation changes when the total number of electrons is kept fixed while changing the 
magnetic field across the saturation value. One can work out easily for free electrons that in order to keep the total 
number of electrons fixed while changing the magnetic field across the saturation value { h = hgat ~^ h^at ~ A/i) one 
has to adjust chemical potential (serving as Lagrange multiplier) accordingly: 

1 TTTI 

Af,{h) = -(4(l-cos— )-A/i) 

(im 4sin2^-A/i\ ^ ^ 

+ 2 cos I — + arcsm ^-J-^ 1 (B4) 

One observes that for all fillings, except for half filling (note that at half filling A^{h) ~ 0): 

^^l{h). 1 



\h=h. 



dh 2 

Calculating the second derivative from (|B4[) we obtain: 

^'^^l{h) _ 1 



dh' l''-^- 2sin2 7rn 



(B5) 



(B6) 
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From Hamiltonian ([3]) we get for the down- spin electrons density: 

< Ani >= -Am = --^/\Afi - Ah/2\ 



(B7) 



Using Eqs. (jB5P and (|B6p we obtain: 



Thus recovering the anticipated result. 



dh 



< An\ >= 



2 9/i2 



Ah 

2tt sin Trn 



= -At 



(B8) 



APPENDIX C: EXACT EXPRESSION FOR SATURATION SUSCEPTIBILITIES 



Here for completeness we write out the Bethe Ansatz results of the susceptibilities for repulsiv o^^i^^ and 
attractivoi^ii^ Hubbard interactions at the saturation magnetization (with the correcting the misprint for attrac- 
tive case). Close to the saturation magnetization ground state energy of the Hubbard model, being an analytic 
function can be expanded in the Taylor series in infinitesimal deviation of the magnetization from its saturation value 
Am: 



where saturation field is 



n 



h , — 9eo\ 
I'' sat ' a™ 



Am) = epM - hsatAi 



(Amf + 



{CI) 



T5^im=n/2 = IhriAm^o --^-^^^^^ — Am ^^"^ invcrsc saturation susceptibility is: 

From the above expansion it follows, that in case ^ 0, A7T1 ^ Ah and in case x^^ = ^ 



\sat 



ra—n/2 • 

Am ~ V Ah. The most general form allowed for the systems where ability to expand ground state energy into the 
Taylor series in magnetization can be proved is Am ~ (A/i)^/", with n being positive integer. 

Repeating the calculation presented in details ii*i^ and correcting the missprint one obtains expressions of saturation 
susceptibilities: 



X-\U<0)\n = -sin(7rn) 



X-HU>0)\n = -sin(7rn) 

TT 



27r — 2 arctan 



4 sin(7rn) 



\U\ 



2 arctan 



4sin(7rn) 
U 



(C2) 



It is easy to verify that these expressions to first order of U indeed recover the bosonization result given by Eq. ^ 
We note as well that in the first order of U susceptibilities for both signs of U have identical form, in accordance with 
the bosonization, this way the typo was identified in the case of attractive interactio n^^i^*^ . 
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